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Abstract 

Electron spin relaxation caused by the D'yakonov-Perel' mechanism is inves- 
tigated theoretically in asymmetrical A3B5 heterostructures. The total spin 
relaxation anisotropy is demonstrated for a wide range of structure parame- 
ters and temperatures. The spin relaxation rates dependences are derived for 
GaAs-based hetero junction and triangular quantum well. The calculations 
show a few orders of magnitude difference in spin relaxation times. 

I. INTRODUCTION 

The degrees of freedom of spin have a great deal of attention throughout the development 
of semiconductor physics. Recently the spin properties of carriers have been investigated 
intensely in low-dimensional semiconductor structures. In electronics, much interest in spin 
has also appeared due to recent proposals to construct spin transistors and spin computers 
based on heterostructures.EH 
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Spin-orbit interaction, which determines a spin behavior, is much more complex in semi- 
conductor heterostructures than in bulk systems. The bulk spin-orbit terms take a more 
interesting form in two-dimensional systems, and, in addition, new terms appear which are 
absent in bulk. 

In Ref. ^ we considered electron spin dynamics in asymmetrical heterostructures. Giant 
anisotropy of spin relaxation times caused by the interference of different spin-orbit terms 
has been revealed. In this work, we calculate spin relaxation rates in real asymmetrical 
structures. A hetero junction and a triangular quantum well are considered in detail. The 
effect of heteropotential asymmetry on spin relaxation is investigated in a wide range of 
electron concentrations and temperatures. We show that giant spin relaxation anisotropy 
governed by external parameters, opening possibility for applications in spin engineering. 



II. THEORY 

Consider a system with spin-orbit interaction described by the Hamiltonian Hso(k), 
where k is a wavevector. Hso(k) is equivalent to a Zeeman term with an effective magnetic 
field dependent on k. In the presence of scattering, the wavevector changes and, hence, 
the effective magnetic field changes too. Therefore for frequent scattering, the electrons are 
subjected to a chaotically changing magnetic field. The spin dynamics in such a system has 
a diffusion character, which leads to loss of the specific spin orientation. This is called the 
D'yakonov-Perel' spin relaxation mechanism,! and it is the main spin relaxation mechanism 
in many A3B5 bulk semiconductors and heterostructures. 

For a two-dimensional system with any Hso(k) (where k lies in the plane of the het- 
erostructure), one can show similarly to Ref. || (see also Refs. f|, |5], |), that the spin dynamics 

of electrons in the presence of elastic scattering is described by the following equations 

00 
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Note that this is true only for times longer than the momentum relaxation time but shorter 
than the spin relaxation times. In (|1|), Si are the spin density components (i = x,y,z), 
the integration is performed over the energy e = h 2 k 2 /2m, where m is the electron effective 
mass, F±(e) are distribution functions of electrons with the spin projection equal to ±1/2, 
Oi are the Pauli matrices. H n are the harmonics of the spin-orbit Hamiltonian: 

H n = <j> H so (k) exp(-in^ k ) , (2) 

where ip-^ is the angular coordinate of k, and the scattering times are 

— = { deW(£,6)(l-coan9), (3) 

where W(e, 6) is the probability of elastic scattering by an angle 6 for an electron with 
energy e. 

Equation ([!]) is valid for two-dimensional electrons with any spin-orbit interaction 
Hso{k). Now we consider an asymmetrical zinc-blende heterostructure. There are two 
contributions to Hso(k). The first, the so-called bulk inversion asymmetry (BIA) term, is 
due to lack of inversion symmetry in the bulk material from which the heterostructure is 
made. To calculate this term, one has to average the corresponding bulk expression over the 
size-quantized motion.! We investigate the heterostructure with the growth direction [001] 
coinciding with the z-axis and assume that only the first electron subband is populated. 
The BIA-term has the form 

H mA (k) = 7 [a x k x (K ~ (k 2 z )) + vyk y ((k 2 z ) - hi)] , (4) 

where we choose x- and y-directions to be along the principal axes in the plane of the 
heterostructure. Here (k 2 ) is the square of the operator (—id/dz) averaged over the ground 
state, and 7 is the bulk spin-orbit interaction constant. It is seen that i^BiA contains terms 
both linear and cubic in k. 

In asymmetrical heterostructures, there is an additional contribution to the spin-orbit 
Hamiltonian which is absent in the bulk. It is caused by structure inversion asymmetry 
(SIA) and can be written asS'iH 
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HsiA(k) = a (a x k y - a y k x ) , (5) 

where a is proportional to the average electric field, E, acting on an electron: 

a = aoeE . (6) 

Here e is the elementary charge and «o is a second spin-orbit constant determined by both 
bulk spin-orbit interaction parameters and properties of heterointerfaces. 

-f^siA also contains terms linear in k. From the (U) it follows that the harmonics with the 
same n are coupled in the spin dynamics equations. This leads to the interference of linear 
in wavevector BIA- and SIA-terms in spin relaxation.! 

For H so = i^BiA + ^siA, the system has the C 2l ,-symmetry. Therefore equations ([!]) may 
be rewritten as follows: 

S z = , S x i Sy = . (7) 

The times t z , t + and r_ are the relaxation times of the spin parallel to the axes [001], [110] 
and [110], respectively. 

If both spin subsystems come to equilibrium before the start of the spin relaxation, then 

F±(e) =F ( f ji ± -e), (8) 

where F is the Fermi-Dirac distribution function and /i± are the chemical potentials of 
electron spin subsystems. If the spin splitting is small, i.e. 

\fji + - a*— I < Ia*+I» l^-l 

then the expressions for the spin relaxation rates 1/Tj {% — z, +, — ) will have the form 





where 
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Eqs. (P|,|T(1D are valid for any electron energy distribution. If the electron gas is degenerate, 
then the spin relaxation times are given by 



— = Ti(k F ) , 



(11) 



where k F is the Fermi wavevector determined by the total two-dimensional electron concen- 
tration, N: 



k F = V2nN . 



(12) 



In this case, the scattering time T\ in Eqs. ( |10"D coincides with the transport relaxation time, 
r tr , which can be determined from the electron mobility. 

For non-degenerate electrons, the spin relaxation times are determined, in particular, by 
the energy dependences of scattering times T\ and r 3 . If T\, r 3 ~ e u , then t^/t\ = const and 
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Here T is the electron temperature and fee is the Boltzmann constant. In the particular 
case of short-range scattering, v — 0, and Ti = r 3 are equal to r tr , which does not depend 
on temperature. 

Spin relaxation times are very sensitive to the relationship between two spin-orbit inter- 
action strengths, ^(k^) and a. From Eqs. flTT|), ( |T3"D it follows that at low concentration or 
temperature, l/r z , l/r_ and l/r+ are determined by the sum of the squares, by the square 



of the sum, and by the square of the difference of j(k 2 ) and a, respectively. This can lead 
to a large difference in the three times, i.e. to a total spin relaxation anisotropy, if 7(& 2 .) 
and a are close in magnitude. 

In real A 3 B 5 systems, the relationship between i^BiA and H$ia may be different. Hbia 
or -£/sia m ay be dominant ,00or they may be comparable .0 

The value of (k 2 z ) depends on the shape of the heteropotential and will be calculated for 
given asymmetrical heterostructures below. The constant 7 is known for GaAs from optical 
orientation experiments.! The correct theoretical expressions for 7 and «o were derived using 
the three-band k ■ p model.B The heterointerfaces give a contribution to cto in addition to 
that of the bulk.El At large wavevectors, a® starts to depend on /cHlI Here we assume 
concentrations and temperatures sufficiently low allowing us to ignore this effect. 

The spin relaxation rates for two types of asymmetrical structures — for a heterojunction 
and a triangular quantum well — are calculated below. The scattering assumed to be short- 
range [y = 0, r 3 = 7~i = r tr ). All parameters are chosen to be correspond to GaAs/AlAs 
heterostructure: 7 = 27 eV-A 3 , m = 0.067m , where m is the free electron mass, «o = 
5.33 A 2 . The time r tr is taken equal to 0.1 ps and assumed to be independent of electron 
concentration. 



III. SPIN RELAXATION IN A HETEROJUNCTION 

In a heterojunction, the size of spin-orbit interaction is governed by the two-dimensional 
carrier concentration, N. (k 2 ) may be estimated as follows:@ 

( k l) = ~ [ —2 J - ( 14 ) 



where k is the dielectric constant. The mean electric field acting on an electron may be 
taken to be equal to half of the maximum field in the junction: 

B=™5. (15) 
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Fig. H] shows the concentration dependence of the reciprocal spin relaxation times for de- 
generate electrons in GaAs/AlAs hetero junction (k = 12.55). The inset shows the spin-orbit 
interaction strengths, ^(k 2 ) and a, and the absolute value of their difference, as functions 
of electron concentration. 

One can see the total spin relaxation anisotropy over a wide range of concentrations. At 
small N, 1/t+ is less than l/r_, whereas at large concentration l/r + is greater than l/r_. 
This is due to the factor which is multiplied by k 2 in fllQ|) being larger for l/r_ than for 
1/t + and the reverse being true for the factors which are multiplied by k 4 . Therefore at a 
certain concentration, the times r + and r_ will be equal. From Eqs. (|10D, (|TTD follows that 
this takes place when 

4 = A{kl) , (16) 

which is fulfilled at AT = 1.1 x 10 13 cm -2 as illustrated in Fig. [|. At larger concentrations, 
spin relaxation is again totally anisotropic. 

Despite the values 7 (A; 2 ) and a are close in magnitude over a wide range of concentrations 
(see the inset in Fig. |I|), all three spin relaxation rates depend on N monotonically. This 
happens because, as the concentration increases, k-p increases as well, and the terms in 
Hso which are cubic in the wavevector become important. The growth of these terms 
with iV dominates the change of the function (a — ^(k 2 )) 2 in (0), hence the monotonous 
concentration dependence of l/r + occurs. 

The situation changes in the case of Boltzmann gas. For non-degenerate electrons, the 
mean value of wavevector and the concentration are independent. For temperatures up to 
300 K, the characteristic k 2 ~ Imk-ol '/h 2 is much less than (k 2 ), and the spin relaxation 
rates are determined by the first terms in (|T3|). As a result, all three spin relaxation times are 
different up to 300 K at a given concentration. The corresponding calculations are presented 
in Fig. g. 

The times r + and r_ are equal to each other at a certain temperature only. According 
to CPL3|) , the corresponding condition is 

7 



Using the GaAs-parameters and v = in ( |14|) , it can be seen that (|17]) is satisfied for 
T « 100 K for AT = 10 11 cm " 2 and at T « 290 K for iV = 5 x 10 11 cm ~ 2 , which agrees 
very well with Fig. ^|. 

At a fixed temperature, the spin relaxation rates are governed by the electron concentra- 
tion. According to Eqs. (p~3|) , the dependences of 1/Vj oniV are similar to the curves in the 
inset in Fig. |I[ In particular, from Eqs. ([Hj) follows that 1/t z and l/r_ have to be close in 
magnitude, and they both greatly exceed l/r + . In addition, l/r + depends on concentration 
non-monotonically. This is confirmed completely by the results presented in Fig. ^. One can 
see that l/r + <C l/r z 1/t_, and the rate l/r+ has a minimum when plotted as a function 
of concentration. This minimum is at N = 1.4 x 10 13 cm -2 , when the terms in ifgo linear 
in the wavevector cancel each other out. The corresponding condition is 

l(k 2 z ) = a. (18) 

At this concentration, the spin relaxation time r + is very large but remains finite due to the 
terms cubic in k. Therefore the difference in the spin relaxation times is more pronounced at 
low temperature. At high T, the terms in ifso cubic in the wavevector become significant, 
and the minimum in l/r + disappears. However l/r + is still much less than l/r_, i.e. 
huge spin relaxation anisotropy occurs in the plane of the heteroj unction even at room 
temperature. 

IV. SPIN RELAXATION IN A TRIANGULAR QUANTUM WELL 

In this Section, we investigate spin relaxation in the following asymmetrical system. We 
consider a structure with infinitely-high barrier at z < and constant electric field, E, at 
z > 0. 

In the framework of this model, 
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(kl) = a [~ ir ) , (19) 

where 

oo 

/ dx [Ai'{x - f3)f 

a = ^ « 0.78 . (20) 

/ dx [Ai(x - (3)f 
o 

Here (—(3) is the first root of the Airy function: 

Ai(-P) = , p » 2.338 . 

The value of a is given by @. 

In Fig. [| the spin relaxation rates are plotted for the triangular GaAs quantum well at 
different electric fields. It can be seen that the total spin relaxation anisotropy occurs for both 
degenerate and Boltzmann gases in wide ranges of concentrations and temperatures. The 
times t + and r_ coincide only at a specific concentration or temperature. For degenerate 
electrons, according to Eq. (|16|), the corresponding curves have an intersection at iV ~ 
3.4 x 10 12 cm ~ 2 for E = 10 5 V/cm and at iV « 7 x 10 12 cm ~ 2 for E = 3 x 10 5 V/cm 
in agreement with Fig. |4]a. For a Boltzmann gas, the intersection of r + and r_ occurs at 
T w 150 K for £ = 10 4 V/cm and at T w 240 K for E = 2 x 10 4 V/cm, according to (0). 
This is also confirmed by Fig. ^jo. 

The behavior of the reciprocal spin relaxation times in electric field is presented in Fig. |5| 
for both degenerate and Boltzmann electron gas. The dependences of 7(fc 2 ) and a on electric 
field are similar to the inset in Fig. [I]: their values are close to each other in magnitude, so 
the difference between them is very small. This leads to the minimum in the dependence of 
1/t + on E. The cancellation condition flUf ) is fulfilled at E w 1.9 x 10 6 V/cm. The electric 
field of this size can be created in heterostructures containing a gate, allowing experimental 
observation of non-monotonic spin relaxation rate dependence shown in Fig. |5|. 
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V. CONCLUSION 

It has been showiMH that inclusion of both BIA and SIA terms (ffl) and (|5]) into i?so 
leads to the conduction band spin-splitting anisotropy in fc-space in A 3 B 5 semiconductor 
heterojunctions. However the spin relaxation analysis that has been performedEl ignored 
this effect. 

The authors of Ref. ^ showed that the BIA and SIA terms interfere in weak localization 
but are additive in spin relaxation. In this paper, we prove that the terms in ifgo linear in 
the wavevector cancel each other in spin relaxation as well. 

In a recent experiment .0 the spin relaxation anisotropy was observed for non-commonly 
used (110) GaAs quantum wells. In this experiment, the spin relaxation in the growth 
direction was suppressed because of the "built-in" anisotropy of the sample due to presence 
of heterointerfaces. In the present paper, we predict spin relaxation suppression in the plane 
of a heterostructure. Moreover, all three spin relaxation times are different in our case, and 
this effect takes place in ordinary (001) heterostructures. 

To observe the predicted spin relaxation anisotropy, one can perform time-resolved mea- 
surements similar to Ref. In steady-state experiments, spin relaxation can be investigated 
via the Hanle effect. To obtain the spin relaxation times, one has to take into account of 
the fact that the Lande g-factor has not only diagonal in-plane components (g xx ) but also 
off-diagonal ones (g X y) in asymmetrical heterostructures.il The degree of polarization of the 
photoluminescence in a magnetic field B _L z is described by the following expression 

P ^ = i _i_ i ( P }°\ R12 > ( 21 ) 

1 + [[IB {gxx ± Qxy) B\ 2 T z T T 

where the upper and lower signs correspond to the experimental geometry B||[110] and 
B||[110], respectively (/zb is the Bohr magneton). 

We show that the terms in the spin-orbit Hamiltonian which are linear in the wavevector 
interfere which leads to the huge anisotropy of the spin relaxation times. At a high con- 
centration or temperature, this effect starts to disappear due to domination of the terms in 
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Hso cubic in k which are present only in Hbia- However the higher order terms in f/siA 
are not forbidden by symmetry either. These terms can also interfere with these in i^eiA, 
and cause additional non-monotonic peculiarities in spin relaxation times dependences on 
the structure parameters. 

In conclusion, we calculated the spin relaxation times for an A 3 B 5 heteroj unction and 
triangular quantum well. The observance of total spin relaxation anisotropy is predicted in 
a wide range of structure parameters and temperatures. 
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FIGURES 

FIG. 1. The concentration dependences of the reciprocal spin relaxation times, l/r 2 (solid line), 
1/t_ (dashed line) and l/r+ (dotted line), for GaAs/AlAs heterostructure at zero temperature. 
The parameters are given in the text. The inset shows the spin-orbit interaction strengths, j(k%) 
(solid line), a (dashed), and \^{k1) — a\ (dotted), in eV-A, as functions of electron concentration, 
N/(10 12 cm" 2 ). 

FIG. 2. The temperature dependences of spin relaxation rates, l/r_ (solid line), 1/t z (dashed 
line) and l/r+ (dotted line), for GaAs/AlAs heterostructure at different electron concentrations. 

FIG. 3. The concentration dependences of the reciprocal spin relaxation times, l/r+ (solid line), 
l/r 2 (dashed line) and l/r_ (dotted line), for Boltzmann electron gas in GaAs/AlAs heterostructure 
at different temperatures. (1) T = 30 K, (2) T = 77 K, (3) T = 150 K, (4) T = 300 K. 

FIG. 4. The spin relaxation rates, l/r 2 (solid line), l/r_ (dashed line) and 1/t+ (dotted line), 
in the triangular GaAs quantum well at different electric field, a — degenerate electron gas, b - 
Boltzmann gas. 

FIG. 5. The spin relaxation rates, l/r+ (solid line), l/r 2 (dashed line) and l/r_ (dotted line), 
in the triangular GaAs quantum well as functions of the electric field, a — degenerate electrons, 
(1) N = 10 11 cm" 2 , (2) N = 3 x 10 11 cm~ 2 , (3) N = 5 x 10 11 cm" 2 , (4) N = 10 12 cm" 2 , b - 
Boltzmann gas, (1) T = 30 K, (2) T = 77 K, (3) T = 150 K, (4) T = 300 K. 



14 




5 10 15 20 25 



Electron concentration, 10 12 cm 



Fig.1 Averkiev et al. 
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Fig. 2 Averkiev et al. 
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Fig. 3 Averkiev et al. 
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Fig. 4a Averkiev et al. 
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Fig. 4b Averkiev et al. 
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Fig. 5a Averkiev et al. 
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Fig.5b Averkiev et al. 



